Fluctuation relations and rare realizations of transport observables 
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Fluctuation relations establish rigorous identities for the nonequilibrium averages of observables. 
Starting from a general transport master equation with time-dependent rates, we employ the stochas- 
tic path integral approach to study statistical fluctuations around such averages. We show how under 
nonequilibrium conditions, rare realizations of transport observables are crucial and imply massive 
fluctuations that may completely mask such identities. Quantitative estimates for these fluctuations 
are provided. We illustrate our results on the paradigmatic example of a mesoscopic RC circuit. 



In the past decade, the concept of fluctuation rela- 
tions has become a powerful new paradigm in statis- 
tical physics. Generalizing the celebrated fluctuation- 
dissipation theorem, fluctuation relations establish con- 
nections between the nonequilibrium stochastic fluctu- 
ations of a system and its dissipative properties. The 
ensuing theoretical and experimental perspectives have 
sparked a wave of research activity [THS] . 

Fluctuation relations generally relate to observables 
carrying thermodynamic significance, such as work, heat, 
entropy, or currents. By way of example, consider the 
Jarzynski relation 6J 



( 



e ^^) = e-^^^. 



(1) 



where W is the work done on a system during a nonequi- 
librium process in which driving forces act according to 
some prescribed 'protocol'. The averaging (...) is over 
all realizations of the process, AF is the free energy dif- 
ference between equilibrium states with forces fixed at 
the initial and final value, and /3 = T^^ is the inverse 
temperature (fc^ = e = 1 throughout). Relations ('theo- 
rems') of this type are usually proven with an incentive to 
establish rigorous bounds on thermodynamic quantities. 
However, considerably less efforts fT", "8] have been put 
into quantitatively exploring the fluctuations of observ- 
ables around the constraints imposed by fluctuation rela- 
tions. For example, for protocols returning to the initial 
configuration, AF = 0, the Jarzynski relation has the 
status of a sum rule, {X) — 1, for the average of the sta- 
tistical variable X = e~^^. Since the average growth of 
entropy demands (}¥) > 0, the sum rule quantifies the 
presence of exceptional process realizations with M^ < 
[51 |S] . The exponential dependence of the variable X on 
W implies that this variable acts as a 'filter' whose fluc- 
tuations around the unit average contain specific infor- 
mation on rare processes. The statistics of X is also vital 
to the applied relevance of Eq. ([l]): as we demonstrate 
below, even modest departures off thermal equilibrium 
tend to amplify fluctuations to the extent that the aver- 
age {X) = 1 cannot be resolved in practice. Conversely, 
observables like X may serve as diagnostic tool detect- 
ing how far away from equilibrium the system has been 



driven, and what physical processes are responsible. 

In this work, we provide a theoretical approach to the 
rare event statistics of fluctuation relations in transport. 
The motivation for emphasizing transport lies in the ap- 
plied importance of 'particle current flow' to the nonequi- 
librium dynamics of stochastic systems. (To give just two 
examples, let us mention charge currents in electronic 
circuits and the migration of species in biological sys- 
tems.) Conceptually, the current flow through a system 
acts as a source of noise on account of the discrete nature 
of particle exchange. Off equilibrium, this 'shot noise', 
rather than thermal noise, is often the dominant source 
of stochasticity, and a consistent theory must account for 
the ensuing feedback cycle of current into noise and back. 
To this end, we will employ Markovian transport mas- 
ter equations [T0Hl2] as a minimal framework resolving 
the statistics of individual particle transmission events. 
We analyze the master equation in terms of a stochas- 
tic path integral [121 US] , an approach tailor-made to the 
description of rare events in analytic terms. In applica- 
tions to 'mesoscopic' systems, the path integral affords 
an interpretation as the semiclassical limit of a quantum 
nonequilibrium Keldysh theory [T^l [H] , thus establishing 
connections between classical and quantum fluctuations. 
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FIG. 1. (Color online) Left panel: Schematic illustration of a 
generic system exchanging particles with M reservoirs. Right 
panel: Realization in terms of an electric RC circuit. 

Generally speaking, we are interested in the statistical 
properties of particle currents, Iu=i,...,m, flowing from 
some 'system' to M connected 'reservoirs', see Fig.[TJ As- 
suming particle number conservation, the instantaneous 



number of particles in the system n is subject to a conti- 
nuity equation, (n) +X]i/(^>') = *-*• The system exchanges 
particles with the t^th reservoir at time-dependent rates 
g^j, see Fig. fll We require a detailed balance condition 
to hold, 






-P-^An), K,,tin)^dnUin)-U,t, (2) 



stating that the logarithmic ratio of rates is governed 
by a cost function measuring the difference in 'ener- 
gies' E^{n) = U{n) — nf^ before and after a parti- 
cle has entered the system through the i^th terminal, 
T\n{g+/g-) = -[E,{n + 1) - E^in)]. Here, U is a 
time-dependent driving force due to reservoir ly, and U{n) 
measures the system's internal energy, where nonlineari- 
ties in U{n) describe particle interactions and dnU{n) = 
U{n + 1) — U{n). The probability Pt{n) to find n par- 
ticles at time t is then governed by the one-step master 
equation pBfl2] 



dtPt{n) = -HgPt{n), 



(3) 



with 'Hamiltonian' Hg(n,p) = '}2v ±\^ ~ ^^^) dvi.'^)^ 
where e^ (e~^) raises (lowers) n by one unit. At the 
initial time t = — r, the system is assumed to be in equi- 
librium, /^,_^ = and P-r{n) = p{n) = e-'^C^C")-'^) 
with free energy F — — Tin ^„ exp(— /?[/). During the 
time interval [— r, r], the external forces f^^t then evolve 
according to a prescribed cyclic protocol, returning to 
fv,T = at the final time t — t. Equations ^ and (Is]) 
describe a large spectrum of transport processes in and 
outside physics |15] . Examples include charge transport 
in mesoscopic devices [IT], molecular motors [I6j, chem- 
ical reaction networks [17], and evolution in biological 
quasispecies models [T5]. 

The master equation ([3]) deliberately emphasizes the 
analogy to an imaginary-time Schrodinger equation with 
Hamiltonian Hg. This formal correspondence implies 
that the time evolution of Pt{n) can be represented in 
terms of a path integral O [13]. Applying a stan- 
dard Trotter decomposition to the unit-normalized 'par- 
tition function', Z = X]n^'^("') ~ ^^ '^^^ obtains 
Z = J D{n,p) e^-^'^'^P"-"'^"'^'^'^ p{n^r), where the in- 
tegration is over smooth paths {(nt-,pt)} and the 'mo- 
mentum' pt G iM conjugate to the particle number is 
integrated over the imaginary axis. For technical details 
concerning discretization and normalization issues, see 
also Ref . |19j . To extract information on the current pro- 
files {/i/,t}, we introduce sources, Z — > Z[x], by coupling 
Hg to 'vector potentials' (counting fields) x — {Xi^,t}, 



Hg{n,p) -^ Hg{n,p,x) = XI ( 



i^,± 



1-e 



Tip-lX:^ 



9y{n)- 



Much like in ordinary quantum mechanics, moments of 
the currents can then be generated by functional differ- 



entiation, 



i5 



SXl^uti 5Xu2,t2 



x=o 



z[x] 



^Z[x] = {e-'^'^^-^'''^^'^). (4) 

This identifies the sourceful partition function 

Z[x] - J D{n,p) e-^«["^P^>^V("-r), (5) 

Sg[n,p,x\^ - dt{pn- Hg{n,p,x)) 

as generating functional. The current probability distri- 
bution function, P[I] = P[/i, . . . , /jv/], follows by func- 
tional integration. 



P[I] 



Dx e 



iT-.F-^dt X.I. 



z[x]- 



(6) 



Following general principles [2D], we now aim to relate 
the functional Z , describing evolution governed by rates 
g = {(7i/,t}, to the functional Z;, computed for time- 
reversed rates {Tg)t = g-t, i.e., Z^ = Z\ ^^ . Here we 

have defined a time reversal operator T acting on 'scalar' 
functions x = {n,g,f) as {Tx)t — x_j, while 'vectorial' 
functions v — {I,p, x) transform as {Tv)t — —V-t- Using 
Eq. (2), it is straightforward to verify that the action in 
Eq. ( 5 1 exhibits the invariance property 

Sg[n,p,x]^Sfg[fn,f{p~pdnU),f{x + tPf)] + 

+ /3[Uinr)-Uin^r)]. (7) 

Inserting Eq. ([7]) into Eq. ([s]), and changing the integra- 
tion variables {n,p) — ^ {Tn,Tp), we obtain the symme- 
try relation Z[x] = Zi,[T{x + if^f)]- Substituting this into 
Eq. ^ , we arrive at a variant of the Crooks relation [3T] , 



P[I] 
Pb[fl] 



-fii:^dtj:^ui^ 



(8) 



where Pb is the current probability distribution computed 
for time-inverted rates Tg. Equation ([8| was first derived 
in Ref. 20] by considering the symmetry of the operator 
Hg generating the system's Markovian dynamics. Our 
present derivation has the advantage that it is based on 
a path integral representation. This gives us the option 
to explore the fluctuation statistics of currents beyond 
the rigorous bound imposed by Eq. ((sl). To elucidate 
this point, let us integrate Eq. ([8]) over / and use nor- 
malization, J DIPi,[TI] = 1, to obtain a variant of the 
Jarzynski relation 



(X) = l, X = e^^-^'''^-^f-^''^ 



(9) 



On average, the currents (1,^) ^ ~Pgufu follow 
the driving forces, which means that typically X ~ 
g-2r/3 }^,^g.f^ jg g^j^ exponentially small quantity. The 



lllll(4^^^^^^^^^^ 



law distribution. This indicates a divergence of all mo- 
ments of X except for the first, (X'^^) = 1. 




FIG. 2. (Color online) Fluctuations of the variable X, Eq. (l9|, 
in the RC circuit over 1000 simulated process runs. For 
tv/tr = 100 (upper panel), the average {X) = 1 remains 
visible. However, for ty /tx = 1 (lower panel), fluctuations 
dominate and conceal the unit average. 



average value {X) = 1 is due to exceptional processes 
where currents fluctuate against the driving forces [9] . In 
this sense, X acts as a selective observable filtering rare 
events. Note that the exponent in X has a clear physical 
meaning: it refers to the work done on the system during 
the process. 

Before discussing the statistics of the variable X in gen- 
eral terms, let us consider the example of a mesoscopic 
RC circuit biased by an external voltage, /i/2,t = ±Vt/2, 
see Fig. IT] While in earlier theoretical [55] and exper- 
imental [23] studies, fluctuation relations for circuits of 
this type have been discussed for the thermal-noise dom- 
inated regime, we here consider the more general case of 
noise self-generated by transport out of equilibrium. The 
circuit's stochastic evolution is described by the path in- 
tegral ([5]) with M = 2 and sequential tunneling rates 



m, 



5^4 («) 



1 



±Ky^t{n) 



R ei/S^-tW-l' 



n,,t ^ d,,U + {-Y 



Vt 



(10) 

where U(n) = §Q{n — 1/2)^ is the charging energy (as- 
suming Coulomb blockade peak conditions). We have 
performed numerical simulations of Eq. ([3J) to compare 
with the analytical estimates described below. Assum- 
ing Ri ~ R2 ~ R, the dynamics is characterized by 
a number of time scales, including the RC relaxation 
time tuc = RC/2, the inverse of the mean rate at which 
charges enter the central island ty = R/V, the scale at 
which thermal fluctuations lead to system-reservoir par- 
ticle exchange, tx = R/T, and the typical scale t^ of vari- 
ation of the external voltage protocol. For tv/tr ^ 1, 
we are in a thermal regime where noise is of Johnson- 
Nyquist type and fluctuations are comparatively benign. 
However, the opposite shot noise regime, tv/tr ^ Ij is 
governed by much more aggressive fluctuations, as illus- 
trated by the simulation results in Fig. [2] Figure [Sj shows 
results for the ensuing changes in the probability distri- 
bution P{X). Comparing the two profiles, one notices 
the massive broadening of P{X) upon entering the shot- 
noise dominated regime. The near linearity of P{X) on a 
double-logarithmic scale suggests a crossover to a power- 




FIG. 3. (Color online) Main panel: Probability distribu- 
tion of X in Eq. ([9| obtained from 2 x 10® simulation runs 
for the asymmetric voltage pulse Vt = Vot^it/il? + t^), with 
t^ = O.Olr, VoT = 250, t/Irc = 2 x lO". The solid black 
curve (Vo/T = 0.5) and the dashed red curve [Vq/T = 3) 
probe the thermal and the near shot-noise regime, respec- 
tively. Lower inset: Simulation results for {X^) vs Vqt/T for 
U, = O.Olr, VoT = 5, t/Irc = 2 X 10*. The straight black hue 
is a fit to Eq. (Ill with I — 2. Upper inset: The black curve 
gives a typical current profile for t^ = 0.04r, Vqt — 500, 
r/tRc = 2 X 10^ G = i2R)~\ and Vo/T = 10^ The red 
curve denotes the averaged profile. 

To make these observations more quantitative, we next 
consider the moments (X'-^) within the path integral 
approach. While details of the discussion are adjusted 
to the above circuit example with AI — 2 terminals, the 
overall strategy is general and can readily be adapted. 
Equations ^ and ^ imply that (X') = Z[i(3lf]. The in- 
troduction of sources renders the minimum of the action 
Sg[n,p,i(]lf] non- vanishing. For sufficiently large obser- 
vation time 2t, the action is large and stationary-phase 
methods apply. We solve the corresponding Hamilton 
equations, fi = dpHg and p = —dnlig, under the as- 
sumption that the driving forces vary over time scales 
tj larger than the intrinsic relaxation times. Again, we 
start from an equilibrium state and assume a cyclic force 
protocol, f±T = 0. Under these conditions, there ex- 
ist quasistationary solutions ?T.,p ~ 0, and the equations 
become algebraic. The further procedure sensitively de- 
pends on whether we are in the thermal regime defined 
by /3/ < 1 for all times, where /( = |/i — /2I, or in the 
off-equilibrium case realized otherwise. 

In the thermal regime, fluctuations around the av- 
erages (n) and (I^) are moderate, and the path inte- 
gral can be expanded to second order in the variables 
p + 131 f^. This expansion, equivalent to the Kramers- 



Moyal expansion of the master equation [TU' , reduces the 
path integral to the Martin-Siggia-Rose functional cor- 
responding to the Fokker-Planck equation [12]. In this 
regime, we may employ the high-temperature expansion 
of Eq. (jil, 5^ w ^1,(1 =F ^)- While the ensuing equa- 
tions, quadratic in p, afford straightforward solutions, 
details depend on the n-dependence of the rates. We 
here consider the prototypical situation of n-independent 
equilibrium rates g^. It is then straightforward to show 
that p = —W^^^fi^g^f^^g,^, while n is determined 
implicitly by the condition h = — (/i) — {I2) = ^ 
^,j9iyKi, = 0. Substituting this configuration into the 
action, we obtain 



{X^ 



exp 



l{l-l)P dt{I)f 



where (/) 



-(/i 



/? 



5192 

91+92 



(11) 



/ is the average time- 



dependent linear-response current (we assume /i > /2). 

In the nonequilibrium regime, however, the shift x = 
ipif represents a massive (exponential) intrusion into the 
action. The stationary phase equations are then solved 

by the optimal p-configuration p = ^ In ( ^" +ig/„ - ) > 
while the sub-exponential action dependence on n is less 
important. Substitution of {n,p) into the action leads to 
a super-exponential scaling with the driving force [24j, 



(X') 



exp 



dt{I) (- 



£2 



{l~l)f 



(12) 



This result shows that the moments do not diverge but 
become extraordinarily large upon entering the nonequi- 
librium regime. Once the condition /3ft <C 1 is violated, 
the unit average {X) = 1 is completely masked and the 
fluctuation relation ^ looses its practical meaning. 



Equations ( 11 ) and ( 12 ) exemplify how the fluctuation 



statistics of X changes dramatically upon departing from 
equilibrium and how these changes contain telling infor- 
mation on the relevant nonequilibrium processes. For ex- 



ample, Eqs. ( 11 ) and ( 12 ) have been derived for the rates 



(10), which in turn rely on the assumption of a uniform 



temperature T (generally established by external cool- 
ing.) In the complementary case of thermal isolation, 
the system will heat up by the very currents the FRs are 
probing. In this case, the statistics oi X - which should 
be straightforward to measure experimentally - couples 
to the effective particle distributions building up in the 
system. Comparison to analytic results may then probe 
the validity of models for nonequilibrium transport and 
the ensuing theoretical descriptions. It is in this sense 
that we believe the fluctuations of X (and other vari- 
ables entering FRs) to contain far-reaching information 
beyond the rigorous FRs themselves. 
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